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Abstract 

It is shown that the four dimensional antiferromagnetic lattice </> 4 model 
has the usual non-asymptotically free scaling law in the UV regime around 
the chiral symmetrical critical point. The theory describes a scalar and a 
pseudoscalar particle. A continuum effective theory is derived for low energies. 
A possibility of constructing a model with a single chiral boson is mentioned. 

I. INTRODUCTION 

This is the second of two papers where we study the impact of higher derivative terms 
in field theories. In ref. [0 we have pointed out that the presence of these terms in a self 
interacting single component scalar field theory produces tree level effects which may drive 
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the formation of new vacua which is not accounted for by the decoupling theorem || . Three 
phases have been studied in the $ 4 theory by means of the mean-field approximation, the 
paramagnetic, < >= 0, the ferromagnetic, < >= $! ^ and the antiferromag- 
netic where < Q(x) > is an oscillating function. 2 d bands have been found in the dispersion 
relation for the elementary excitations above these vacuua in a certain range of the coupling 
constants of the d-dimensional theory. A reduced Brillouin zone was introduced for each 
band. Two zones describe particle like excitations and the others decouple in the mean-field 
continuum limit. 

For a special choice of the coefficients of the higher order derivatives the theory possesses 
a formal chiral symmetry which allows us to decouple the two particles. The two decoupled 
modes correspond to the sublattices consisting of the even and the odd lattice sites. The 
theory which has nonvanishing field variables on one of the sublattices only is local and 
describes a chiral boson. In fact the space inversion exchanges the two sublattices and there 
is no space inversion partner of the particle in such a model. 

We extend the analysis |1J in the upper critical dimension, d = 4, for theories in the 
vicinity of the chiral invariant critical point to the one-loop order and show that the beta 
functions of the lattice regulated theory with our 0(D 2 ) term in the lagrangian are those 
of an ordinary $ 4 model and give a renormalized lagrangian in terms of the continuum 
field variables. The one-loop renormalizability turns out to be a nontrivial consequence of 
the perturbative renormalizability around the critical point of the usual $ 4 model without 
higher order derivative terms because we have to render the dynamics for the two particles 
finite by fine tuning the set of the parameters of the bare lagrangian with a single quantum 
field. 

There is a formal similarity between the tricritical point of the 6 model and the chiral 
symmetrical theory. The mean field solution of the model with the potential 

n0) = |0 2 + |0 4 + |0 6 (1) 

shows a tricritical point at = which separates the second and the first order phase 
transition lines with different scaling laws ||. In our case the dispersion relation 

G-\p) = m 2 +p 2 - C2 a Y + c 4 a 4 p 6 (2) 

produces a tricritical point when c 2 is sufficiently large to give an absolute minimum at 
nonvanishing values of the momentum. When G~ l at the minimum is negative an inho- 
mogeneous condensate is formed. The chiral symmetrical point where this happens is the 
Lifshitz point. This was introduced in Q where the e-expansion was used to find out the scal- 
ing laws. A scalar model where the dispersion relation has a single minimum at nonvanishing 
momentum was considered in || and 0. The phase transition towards the inhomogeneous 
vacuum was identified in the mean field level and the quantum fluctuations were taken into 
account in ||. We will be working at d = 4 and extend the loop computation into the phase 
with inhomogeneous condensate which generates a "dangerous irrelevant variable" ]/]]. The 
dispersion relation of our model has several minima hence it contains several particle modes 
simultaneously. The condensate formation mechanism selects one of these particle sectors 
in a manner reminiscent of the spontaneous symmetry breaking. 

There is a technical problem to solve in achieving this goal because more than one particle 
corresponds to the same quantum field. The formal problem is that higher order derivative 
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terms in the kinetic energy imply the presence of states with negative metric || and may 
render the effective action complex. But we argue that far below the momentum scale of the 
condensate we find only two particles, both with positive metric. Their dispersion relations 
can be replaced in the continuum limit by the usual quadratic expressions coming from a 
manifestly hermitean free lagrangian ||. The situation turns out to be somehow similar to 
the species doubling of the lattice fermions where one finds several particle modes in the 
dispersion relation of a single bispinor field. We introduce a 2 d -component field variable, 
<& a (x), a = l,---,2 d , for the computation of the one-loop generator functional for the 
1PI functions of the different excitation bands and show that it can be made finite by an 
appropriate fine tuning of the coupling constants of the original lagrangian. 

The organization of the paper is the following. The basic tools of the perturbation 
expansion are developed in Section 2. The computation of the effective potential is presented 
in Section 3. The elimination of the divergences is shown and the finite renormalized coupling 
constants are obtained in Section 4. Section 5 is devoted to a simple effective theory which 
reproduces our model at low energy. A brief conclusion is in Section 6. 



II. THE PERTURBATION EXPANSION 

We develop the basic formula for the perturbation expansion in the scalar c/> 4 model in 
d — 4 with higher order derivatives by keeping the original field variable, As in |l[ the 

theory is regularized on the lattice. 



A. The lagrangian 

The model considered contains a one component field variable, and is defined by 

the bare, cut-off lagrangian, 



C = -^$(x)/C^n)^$(x) + ^$ 2 (x) + ^$ 4 (x), (3) 
where 

K{z) = 1 + c 2 z. (4) 

We write this lagrangian as 

C = d + C 2 (5) 

with 

A = Xfnk t + ^ + ^, (6, 



and 



£2 = — <9„$/C ^— -□ )d u $+— — $ 2 + — $ 4 . 



(7) 
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We will use C\ non-perturbatively in the selection of the saddle point and C 2 perturbatively 
in removing the UV divergences of the loop-corrections. The bare parameters are defined as 
m 2 B = m 2 R + 5m 2 and A# = \r + 5X. We have no counter terms for the coupling constants 
Cj because their leading order renormalization is at the tree-level. 

We employ lattice regularization where one introduces the dimensionless variables x^, 
(p = a d / 2_1 $, m 2 L = m 2 a 2 and the unit vectors (e^)" = 5^ and writes the action as 



where 



Si[f] = £A(x) = Y\k[wA + V(<p(x))V (8) 

X X ^ ' 

K[ip; x] = ~(p(x) A^(x) + ( B (v( x + e n) + <P(x - e /J) 
+ C(<p(x + 2e M ) + <p(x - 2e„)) 

+ (e(ip(x + e M + e v ) + 2ip(x + e M - e v ) + <p(x - e M - (9) 



and 



= ^V(s) + ^V(s), (10) 



2 r v / 4 

The coefficients of the kinetic energy are 

A = -2d + (4rf 2 + 2d)c 2 , 
£ = 1 - 4cfc 2 , 
C = c 2 , 

E = c 2 . (11) 
The field variable, ip = ip vac + 0, is the sum of the tree-level vacuum, 

(f vac (x) = ipi + (f 2 ix(x), (12) 

where 

d 

X(x) = (-1)»= 1 , (13) 

and the quantum fluctuations, <f>(x). We will study the theory in the para- (ipi = (p 2 d = 0), 
ferro- (ipi 7^ 0, ip 2 d = 0) and the (1, 2) antiferromagnetic (<fi = 0, ip 2 d 7^ 0) phases in 
d — 4. The lagrangian for the quantum fluctuations is 

C 1P = I^(x)X:(^n)^(x) + + ^0 4 (x), 

C 1F = I^(x)/C(^n)^(x) + + 3A^ 2 )0 2 (x) 

+A^i0 3 (a:) + ^V(x), 
Aaf = ^ M 0(x)/C(^n)^0(x) + i(m| fl + ?,\ R ^ 2 d)(t> 2 {x) 

+X RV2 d X (x)<j ) 3 (x) + ^V(x), (14) 



C 2P = ^^(x)/c(^n)^(x) + ^\x) + ^V(*), 
£ 2 f = ^^0(x)/c(^n)^0(x) + Vl {8m 2 L + 5A^)0(x) 

Aaf = ^9 M 0(x)/C^^-D^9 /1 0(x) + ip 2 dx(x){6m 2 L + 8\ipl d )<j)(x) 

1 <5A 
-{5ml + 35A^)0 2 (^) + 5Av9 2d x(^)0 3 (^) + ^0 4 (^), 

where the tree-level vacuum is given by 

P : tp P1 = Lp P2 d = 0, 
F : (p F1 = -^(p F2d = 0, 

m\ R + .M| 

: (p AF1 = <^AF2d = ^ • 

Here M 2 L stands for the eigenvalue of the the kinetic energy on the antiferro magnetic 

with 

M 2 L (d, c 2 ) = 4d/C(-4d) = 4d(l - 4dc 2 ). 

The tree-level conditions for the three phases shown in Fig. 1 are 

P: m 2 LR >0 m 2 LR + M 2 L >0, 
F: m 2 LR <0 M 2 L >0, 
AF : m\ R + M\ < 0M 2 L < 0. 



B. The free propagator 

The free propagator, 
is given by 

G-\p) = m 2 LR +p^lC(-p^), 
where mass parameter with the tilde includes the shift due to the condensate 



Z,2 



m 



LR 



-2ml R F, (22) 

{-2ml R -3M 2 L (d,c 2 ) AF, 



in the different phases and 

/;,. = 2 sin 



p„ = 2sin^. (23) 



We further write 



with the help of 



G-\p) =V\ V ) - c 2 V\p) + m\ R) (24) 



f 1 

It is advantageous to divide the Brillouin zone 



7> 2 (p)=4£sin 2 ^-. (25) 



into 2 d restricted zones, 



whose centers are at 



B=\k„, |M<tt}, (26) 



B a = \\k,-P,(a)\<^\, (27) 



P^a) = irn^a), (28) 
where n^(a) = 0, 1 and the index 1 < a < 2 d is given by 

= l + j^a^-\ (29) 



a 



The propagator for the zone B a is 

G a (p)=G(P(a)+p). (30) 

It turns out that only the Brillouin zones a — 1 and 2 d contain particle like excitations and 
the corresponding propagators are 

G-\p) = rh 2 LR {a) + Z(a)p 2 + 0(p 4 ), (31) 

where the mass and the wave function renormalization constant are given in Table 1. Note 
that m 2 LR (l) =m\ R . 

The fact that the vacuum is a single Fourier mode offers the possibility of recovering the 
energy-momentum conservations in the anti-ferromagnetic phase. The possible translations 
which keep the vacuum invariant consist of an even number of shift of the integer lattice 
coordinates. The corresponding spectrum of the momentum operator is 

Pafu = (modvr). (32) 

In fact, the function mod7r substracts the part of the momentum which can be exchanged 
with the antiferromagnetic vacuum and the resulting value is conserved. In this manner 
the momentum non-conservation is traded for the exchange of the particle type, the "flavor 
dynamics" . 
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C. Chiral symmetry 



The chiral transformation 

X ■ <f>(x) — ► x(x)(/)(x), (33) 

which appears as the shift 

p,->p» + P,(2 d ) (34) 
in the Fourier space is a symmetry of the lagrangian when 

c 2 = c 4 = 0. (35) 

The two particle species are degenerate in the chiral invariant theory. 

The operator V± = \{1 ± x) identifies the fields which belong to the even or odd sublat- 
tices, 

V ± <P ± = 4>±. (36) 

The chiral transformation is represented by 

<t>± -> ±0±, (37) 

so the chiral fields 0+ and 0_ decouple in the chiral invariant theory. The inversion of odd 
number of coordinates exchanges the chiral fields. The low energy excitations in B\ and B\q 
correspond to 

0± = + ±0_, (38) 

where the fields cf>± are slowly varying. Thus the low energy excitations of the zones B\ and 
£>i 6 have space inversion parity +1 and —1, respectively. 



III. THE EFFECTIVE POTENTIAL 

The renormalization of the theory will be performed in the para- ferro and the (1, 2) 
antiferromagnetic phase in the one-loop order by making the effective potential cut-off inde- 
pendent. It is easy to verify that this latter is enough, i.e. the wavefunction renormalization 
constant is finite at the one-loop order, 5Z = 0. 



A. A one-loop diagram 

In order to develop the appropriate notation we consider first a simple example, the 
contribution of the second graph of Fig. 2 in the most complicated case, the (1, 2) antifer- 
romagnetic phase, 
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£(*) = \*r<& J p<7T -0^ G ( k + P)°(P + P ( 2 "))- ( 39 ) 

The lattice cut-off of the loop integrals, p < ir, should always be understood as the constraint 
\Pfj,\ < 7T, for // = 1, • • • , d imposed on the torus V unless it is stated otherwise. An integral 
like this can be written in a simpler form by the help of the following matrix notation. The 
loop integration is split into the sum over the 2 d restricted Brillouin zones, 

/ d d pf(p) = W d d pf(P(a)+p), (40) 

Jp<n a=1 Jp<TT/2 

in particular, 

/ d d P G(p) = £ I d d pf(P(a)+p) = W d d P G a (p). (41) 

Jp<n ^ Jp<7T/2 ^ Jp<n/2 

Returning to our one-loop integral ([39]) we find 

E c (*) = E / , T^Gaik + p)G a (p), (42) 

where we have introduced the region complementer to a, 

a = 2 d + l-a. (43) 

To simplify further the latter expression we now promote a to be an internal index 
distinguishing different kind of fluctuations and define the propagator, 

G a , f3 {p)=5 a , f3 G{P(a)+p.), (44) 

which is diagonal in this new internal space. The contribution considered to the self energy 
is then written in matrix notation, 

nk)*=\#R<&f , 7^tr[G(fc + ph*G(ph 2 ' 1 ], (45) 
by the help of the matrix 

7a,/3 = <W,3,2<*+1> (46) 

which describes the change of the type of particle after scattering off the vacuum. 

B. The one-loop effective potential 

Let us denote the usual 1PI functions by T^ n '{pi, • • • ,p n )- The 1PI function for the 
excitations of the type cui, • • • , a n is given as 

T^(P(a 1 )+p 1 ,---,P(a n )+p n ). (47) 
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The generator function for the zero momentum excitations, the effective potential, is defined 

as 



OO 1 

W*) = £zt E ^-i/'^Nr-^K)). (48) 



71 ' 

n=0 «!,-,«„ 



The matrix 7 2 in (15) reflects a modification of the Feynman rules. Whenever a propa- 



gator G a (p) is inserted in a graph it contains the momentum Pfi(a) + p. We keep track of 
the first term of this sum by introducing a 2 d -component field, $ a , in such a manner that 
the a-th component will be responsible of the excitations in B a . Thus the Feynman rules 
are those of a 2 d -component field with the propagator (f44|) and each external line with p = 
is represented by the insertion of the matrix 

<p=Y,r® a , (49) 

a=l 

where 

d 

Ip,* = IT ^^+a li -p ll (rnodZ),0 (50) 
/x=l 

takes care of the change of the particle type at each vertex due to the momentum flowing 
from the external leg. We will use either the index a or its vector representative, n At (a), in 
the formulae. 

Taking advantage of the matrix formalism introduced above we obtain 

Ve/ f m = U , ^- d tTln[P 2 (P+p)}C(-V 2 (P + p)) 
+ml R + 6X R f^+3X R f), 

= \i , ^rHV 2 (P + P )IC(-V 2 (P + p)) 
2 J p <tt/2 {2ir) a 

+mL + 3A^+0) 2 ], (51) 
where the matrix P is given by 

P a ,p = 5 a ,pP{a) (52) 

and the vacuum field is 

4> = Vil 1 + W7 2d - (53) 
The complete one-loop effective potential V^($) + V^|($) for the background field 

^=$i7 1 + ^7 2d , (54) 

is obtained in (|A3|) and ( |A11|) , 
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and 



^)($) = I^ 2 (P(l))/C(-P 2 (P(l)))+mL + 5mi)^ 



+-(V 2 (P(2 d ))/C(-V 2 (P(2 a ))) + m{ R + 5mi 
+ ^($4 + $ 4 d+6$ 2 $ 2 d)) 



V AF(0) {<5>) 



V p (°\$ + <p AF ), 



Veff ^ 2 W (2tt)^ ^ m l 



2 y p <7r/2 (27T) d 

x V 2 (P(a)+p)IC(-V 2 (P(a)+p)) 
+m 2 LR + 3A«($^ d + $ 2 ) 
V 2 (P(a)+p)K,(-V 2 (P(a)+p)) 



eff 
r AF(l 
'eff 



+m 2 LR + 3A R ($ 2 d 



- 3Q\ 2 R <S> 2 <S> 2 2d 



(55) 



(56) 



The mass parameter of the effective potential in the ferro- and the antiferromagnetic phase 
after the shift $ — > & + cp is given by (^2|) . 



IV. THE RENORMALIZATION IN D = 4 

The divergences arising in the one-loop integral for the effective potential are isolated 
by expanding the logarithm in the integrand. We reintroduce the lattice spacing and use 
dimensional quantities in the rest of this paper. One finds three divergent integrals, 

16 - d^p 



16 „ Jd 

fij= S4(l G " WG,w ' (57) 

and the divergent part of V AF ^> ($) turns out to be 
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where 



and 



AF(1) 



div 



r r 2 



P 9 



-A ii P 1 [($ 16 + ^ 16 ) 2 + ^] 



(58) 



-^A^36P 2 [($ 16 + <^ 16 ) 4 + $ 4 + 6($ 16 + <^i 6 ) 2 $ 21 
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\ 2 R AV, 



($i 6 + <^i 6 ) 4 + $ 4 - -($i 6 + <^i 6 ) 2 $ 2 



AX? 2 = £> 2 - 2£> 2 , 



(59) 



B = 6A R $i(0i 6 + $i 6 ), 

C = 3A fi (2< / ? 16 $ 16 + $ 2 6 + $ 2 ). 



(60) 



The corresponding expression for the ferromagnetic phase can be obtained by the exchange 
1 <-> 16 of the internal index. The condensate has to be set to zero, <p a = 0, in the expressions 
of the paramagnetic phase. 

The choice of the mass counterterm, 



5m 2 L 



-3X R V 1 - 9X 2 rV 2 w V 2 



(61) 



is straightforward after comparing ( 58|) with (|55|) . But there is a problem with the coun- 
terterm <5A because it can not eliminate the divergences for both particles in the same time 
when AX> 2 ^ 0. 



A. Renormalization with chiral symmetry 

The remedy of the problem of the divergences 0($ 4 ) comes from the observation that 
the chiral symmetry protects against the unwanted divergences. In fact, the chiral transfor- 
mation, (0), acts as 

G a (p) - G a (p) (62) 

on the variables of the effective potential and the propagator and the chiral symmetry 
requires 

G a {p) = G & (p) (63) 
which reduces the number of divergences since we gain the relation 

AV 2 = V 2 - 2V 2 = 0. (64) 
The divergent part of the effective potential is now written as 
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^t F(1) ( < f) = 3A i? P 1 [($ 16 + ^i6) 2 + $?] 



- 9A^D 2 [($ 16 + y? 16 ) 4 + n + 6($i 6 + ^i 6 ) 2 ^]. 



Comparing it with (|55|) we arrive at the choice 

SX = l8X R f> 2 . 



(65) 



(66) 



Thus one can eliminate the divergences of the chiral symmetrical theory in either of 
the phases by the help of the appropriate fine tuning of the parameters m 2 B and X R of the 
original lagrangian. The chiral invariant theory is invariant under the exchange of the two 
degenerate particles. Using the chiral fields, $± = $1 ± $i6, one can decouple the two 
particle modes. Let's consider for example the case of the paramagnetic phase. Replacing 
in (^) the appropriate values of 5m 2 L and SX, as given respectively in (|6ll) and (^6|), from 
(|55|) and ( |56"D we get for the effective potential along the chiral line xp ( see Fig- 1)5 



V p „ 

eff 



($1,$ 



16/ 



V C f f 
eff 



($. 



V c f f 

eff 



(<&_ 



(67) 



where 



V^m = lrnl^ + \x R ^ + l - 2 t J <jl 



d 4 p 



x In 



£i J p<^ (2^) 4 
" . • PM ^rw) ■ /') 2 



-3X R V 2 I 
+9X1® 



A 2 
1 



ml + QX R ¥' 



d A p 

P <jL (27r) 4 p 2 + m 2 



(68) 



R 



2 ^ 2 / 1 

' P <£ (2tt)4 (p 2 + m|) 2 - 



The same is true along the chiral lines in the other phases. This decoupling arises because 
in either of the phases at the chiral line the lattice decouples into two different sublattices 
as explained in Ref. ffl. 



B. Renormalization around the symmetrical point 



In order to remove the symmetry with respect to the exchange of the two particles we 
consider the four dimensional theory with the tree-level cut-off dependence 



m 



LR 



2 2 



RU ' ° 2 ~ 16 



1 + a(a/i) 



2 + K 



(69) 



where a = ±1 and \i is a mass parameter to characterize the split of the degeneracy in the 
spectrum, 



M 2 L = -16a(afi) 



2 + K 



(70) 



The quantities referring to the symmetrical theory, fi = 0, will be labelled with a star. Since 
there are several possibilities in reaching the continuum limits as shown in Fig. lb we collect 
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the corresponding conditions for (|69|) in Table 2. We will find that the one-loop corrections 
do not change qualitatively the tree-level spectrum. The degeneracy, m R (2 d ) = m 2 R (l), is 
achieved analytically at the chiral line indicating that the chiral symmetry is not broken 
dynamically. We find m R (l) < m 2 R {2 d ) in the phase F and in the region Pp (Pf is the 
region of the paramagnetic phase on the left side of the chiral line \Pi see Fig. lb). On the 
contrary, m R (2 d ) < m R (l) in the regions AF and Paf (on the right side of the chiral line 
Xp in Fig. lb). We found no singularities in the effective potential due to the discontinuity 
in the momentum of the condensate when the chiral line is reached from the phases F or 
AF. 

The complication we face is that there will be finite /x-dependent corrections from the 
counterterms in the vicinity of the symmetrical theory, T> 2 7^ 2T> 2 . The detailed study of 
the ii dependence in the limit a — > is presented in Appendix fB|. One finds that the in- 
dependence drops out from the finite part of the effective potential and is finite for T> 2 and 
T> 2 . By introducing 

V 2 = T>2 + 5V 2 , 

V 2 = V* + 5V 2 , (71) 
with T> 2 = 2T> 2 one finds the finite expressions 

5V - 1 ln ^ (16) 

^ 2 -"l6^ ln ^(iT' 

1 ^p(l) , ^1(16) 

SV 2 = , k k ' - - In R > / ' 72 

167r 2 m|(16) -m|(l) m|(l) V ; 



with 



and 



m\ P 
m\{\) = { -2ml F (73) 

-2m| + 48/i 2 (a/i) K AF, 

!m R — 16aLi 2 (aii) K P, 
-2m\ + 16 / u 2 (a/i) K F, (74) 
-2m\ + 32LL 2 (aLi) K AF. 

These expressions lead to the counterterms 

Sm 2 = -3X R V 1 - 9X 2 Rl f 2 l6 V 2 , 
5X = 18A|P*. (75) 

It is well known that the spontaneous symmetry breaking in a ferromagnetic theory 
changes the counterterms by a cut-off independent finite piece and influences the renor- 
malization group flow at finite energies only. One could, in principle, encounter a different 
situation in the antiferromagnetic phase because the condensate is formed at the cut-off 
scale. Furthermore one band of the elementary excitations, in B\e, belongs to the staggered 
modes which show fast oscillation at the cut-off scale. It is the fine tuning of the value 
of the minimum of the dispersion relation in the zone Biq which eliminates the divergent 
phase dependence in the counterterms and restricts the effects of the phase transitions in 
the infrared region. 



13 



C. Mass spectrum 



We are now in the position to follow the renormalization in the vicinity of the critical 
system. The effective potential is written as the sum of the finite and divergent part, 



(i) 



eff 



+ V 



(1) 
div i 



(76) 



where the second term in the right hand side is defined by (|58|). One should bear in mind 
that starting with a single mass parameter in the bare lagrangian we have already introduced 
different masses for the propagators in the zones B\ and B\§. The physical masses which 
contain the radiative corrections are given by the derivative of the effective potential. 
The Brillouin zone B\. The mass square of the excitations is given by 



16 



m R + 5m + 3(\ R + 5\)ip 
+3\ R V 1 - 9\ 2 R <p 2 16 V 2 



By the help of the counterterms 



we find 



m ph\ l ) 



m\ + 3A R ^ 6 (1 - ?>\ R 8V 2 - 12\ R 5V 2 ) + dlj^i®) 
—m\ + 16/z 2 (a/z) K 



where 



5V 2 
SV 2 



In 

16tt 2 -m R + 2A^ 2 (anY' 



16tt 2 



-m\ + 16/i 2 (a/i) K j — m 2 R + 16/i 2 (a/i) K 



-8a/i 2 (a/x) K -m\ + 24/i 2 (a/z) K 



(77) 



(78) 



(79) 



(80) 



The computation of the finite part of the effective potential in Appendix C yields vanishing 
result for the second derivatives with respect either field variables in all phases. So we arrive 
at 



■m ph (l) 



in the continuum limit with 



and 



A* 2 («) 



AV 2 (k) 



-2m\ + 48// 2 (/c) - 18A|^ 2 6 A£> 2 (fc), 



lim /i 2 (a/i) K 



lim[P 2 (/i) - 2VM] 

a— >0 




(81) 



(82) 



hm\5V 2 -25V 2 ] 



(83) 



m*(l)-m£(16) 





ln^ 



m|(16) 



K 



0. 

K > 0. 
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The Brillouin zone B\q: One finds 



d ke V e A ff ml =0 = G£(0) + 5m 2 + 3*Atf 6 + 3A fl 2>! 

-36A^ 2 + ^ 16 ^: (1) ($)|^ o 

= -2m\ + 32fi 2 (afi) K + 27\ 2 R ip 2 w (26V 2 - 5V 2 ) 

= m 2 ph (lQ), (84) 

which results 

m 2 ph (16) = -2m\ + 32// 2 (k) - 27A^? 6 AP 2 («). (85) 

Our conclusion is that for k > 8V 2 = 5V 2 = so the two particles become degenerate 
and the chiral symmetry is restored in the continuum limit. For k — the mass spectrum 
stays non-degenerate. Finally the masses diverge as expected when k < 0. 

The ferromagnetic phase: We have, in a similar manner 

m 2 ph (l) = -2m 2 R - 18A^?AP 2 («) (86) 
m 2 ph (lQ) = -2m\ + 16/i 2 ( K ) - 27A^?AP 2 («). 

The paramagnetic phase: The renormalized masses for the line Xp are 

"»i(16) = m% (87) 

showing the presence of the chiral symmetry. In the remaining part of the paramagnetic 
phase we find a non-degenerate spectrum, 

m lh( l ) = m w 

m 2 ph {\Q)=m 2 R + l^ 2 {K). (88) 
The P — AF transition line corresponds to the spectrum 

m lhiX) = m R, 

m 2 ph (16) = 16fi 2 ]imUar. (89) 



D. Coupling constant renormalization 

The Brillouin zone B\ : The definition of the renormalized coupling constant is 



2<T> i & T /AF(1), 



$=0 



$=0 



= 6A ph (l), (90) 

giving 
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1 



With our choice of the counterterms we have 



6 



$=0 



a— >0 r> >fl * J 



1 

6 ^"0 



$=0 



in the continuum limit. 

The Brillouin zone £>i 6 : The self-coupling constant for the field $16 is 



$=0 



V(!6) = A fl - 9A|lim5P 2 + - lim^VjJJW^) 

For the coupling constant which mixes the two fields we have 

d k d L V e A ff^l =0 = 6 ( A * + 5\) - 18A 2 K (P* + SV 2 ) 

-72A^(P* + 8V 2 ) + 91 
= 6A pfc (l,16). 

In the continuum limit it is 

X ph (l, 16) = X R - 3X 2 R \im(5V 2 - A5V 2 ) + \ lim^fl^V^fW^) 

a— >0 h a— >0 



$=0 



$=0 



(91) 



(92) 



(93) 



(94) 



(95) 



The finite part of the effective potential, V^f (<&), is computed in Appendix y. The 
corresponding expressions in the ferro- and the paramagnetic phases are formally the same. 



V. A LOW ENERGY EFFECTIVE THEORY 

Our theory with a single quantum field contains two particles and its antiferromagnetic 
vacuum is in the ultraviolet regime. So it is not obvious that the evolution of the coupling 
constants for the two particle like excitations obeys the renormalization group equations 
which hold for the usual para- or ferromagnetic theories. In order to obtain the renormal- 
ization group equation for the potential of the model we introduce the running cut-off, k, 
implemented in each restricted Brillouin zone in a spherical symmetric manner, 



V a {k) 



<k u 



(96) 



where A = 2n/a and the contributions coming from the edges of the toroidal Brillouin 
zones are left out. We approximate the dispersion relation in T>i(k) and Pie (A;) with an 
0(4) invariant parabola and neglect the non-particle like excitations. These approximations 
involve irrelevant operators of the perturbative continuum limit which should not influence 



the finite energy behavior. Thus the renormalization group equation [ 10 1 for the potential is 
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kd k V£(Qi,$2i) = -a,. 



-" ,k d ln| yk 2 + m|(l) + ^VfW 
k 2 + mU2 d )+dl d V k p (®) 

V k F (<S>) = V k p ($ + ip F ), 



v k AF m = va^+vA F ), 



(97) 



in the leading order of the gradient expansion where stands for the solid angle in d 
dimensions. The coefficients of the higher order terms in k 2 in the logarithm are kept fixed 
in our approximation. 

Consider now the following renormalizable continuum lagrangian for a scalar and a pseu- 
doscalar field, <fi+(x), 4>-{x), respectively, with momentum space cut-off, 



whose renormalization group equation in the leading order of the gradient expansion is 

~ > 1 



(98) 



kduVi 



Q d k d \n\[k 2 + dlV k (4>)}[k 2 + dl V k (4>)} 



(99) 



The renormalization group flow of this model agrees with our higher derivative theory at 
low energies when the initial condition 



ni 



+ 



m 



B— 72 



+ 



(.2 ±2 



(100) 



is chosen. In other words, the model ( |10U| ) is equivalent with (^) at low energy when 
the continuum limit is taken. The correspondence between the phases is the following, 



P >= 0, < 0_ >= 0, 

F 0+ >^ 0, < 0_ >= 0, 

AF 4> + >= 0, < 0_ >^ 0. 



101) 



The conserved momentum of the antiferromagnetic phase is (|32|) and the exchange of 
the momentum 7m(16) M /2 on the lattice with the vacuum corresponds to the exchange 
of the scalar and the pseudoscalar particle. Due to the vertex <\>\§ 2 _ in the lagrangian a 
pseudoscalar particle can decay into two scalar ones in the antiferromagnetic phase and the 
parity is not conserved. 

The one-loop scaling laws of our theory agree with a usual two component <ft 4 model up 
to irrelevant terms. Thus one may suppose that our theory is not asymptotically free and 
consequently becomes trivial in the continuum limit. In this case when the cut-off can not 
be eliminated from the interacting theory the irrelevant terms which were neglected in the 
comparison might be important and generate different physical content. 
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VI. CONCLUSIONS 



The one-loop vacuum polarization effects were studied in the para-, ferro- and (1, 2) 
antiferromagnetic phases of the four dimensional <j) 4 model around the chiral invariant critical 
point. One can identify two particle like excitations in each phase. The one-loop divergences 
were eliminated by an appropriate fine tuning of the parameters of the bare lagrangian and 
the resulting theory was found to be equivalent at low energies with a usual renormalizable 
model made by a scalar and a pseudoscalar field. In this continuum limit where the length 
scale of antiferromagnetic vacuum or the pseudoscalar staggered particle mode tends to zero 
the well known problems about the unitarity || disappear. 

One should emphasize that even though the cut-off can be removed and the continuum 
limit can be taken at the one-loop level the theory can be defined by relying heavily on the 
regulator. The renormalized continuum theory exists only when the regulator is taken into 
account both at the tree- and loop-levels in a systematical manner. 

The possibility of removing the divergences in the presence of an apparently non- 
renormalizable term in the lagrangian is in principle a serious threat for the usual strategy 
of Particle Physics where the universality is used to limit our investigations to the class of 
renormalizable theories. But the result that our model reproduces the infrared structure 
of a conventional renormalizable one is reassuring because it indicates that there is no new 
universality class encountered. 

The antiferromagnetic phase is certainly different compared to the usual <j) 4 model with 
Cj = 0. But even the para- and the ferromagnetic phases in our higher derivative model 
become unusual, as well, in the vicinity of the chiral invariant critical point. This is because 
the dispersion relation develops in all of these phases a second minimum which can be fine 
tuned around this critical point in such a manner that another particle like excitation, the 
analogue of the rotons of superfluids, appears. This particle has staggered excitation modes 
which allow us to introduce the chiral fields which are exchanged between each other under 
the space-time inversions. These chiral fields decouple in the chiral symmetrical theory. 
By considering field variables only in the sublattice of the even lattice one can construct 
models with a single chiral boson. Such a decoupling of the modes is reminiscent of the 
fermion doubling on the lattice and the resulting model with a single chiral boson is local 
and satisfies the reflection positivity. 

Our computation was made at the one-loop level only. There is no conceptual problem in 
extending our work to higher loop orders though the treatment of the overlapping divergences 
with unconventional dispersion relation represents a challenging problem. It remains to be 
seen if the perturbative elimination of the divergences can be achieved beyond the one-loop 
order. If a theory in the antiferromagnetic phase turns out to be renormalizable then its 
vacuum appears homogeneous in physical measurements. It is the structure of the excitations 
only which betrays the non-trivial structure of the vacuum of such a theory. There are 
numerical indications of the continuum limit in the antiferromagnetic phase for other models 
with antiferromagnetic vacuum [ITT . 



One should mention that there are other possible continuum limits in our model away 
from the chiral invariant critical region when the mass parameter is kept at a cut-off in- 
dependent value. This parameter plays a role analogous of the ^-parameter of the Wilson 
fermions. In fact, the excitations of the restricted zone B\ decouple when m 2 LR = O(a ) ^ 0. 
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The only left over excitations in become critical along the P-AF transition line. Thus 
the approach of the critical line with a fixed mass parameter results in a theory which con- 
tains a single pseudoscalar particle. In a similar manner certain regions of the c 4 7^ part 
of the phase diagram may become critical and offer a continuum limit. This is because the 
renormalizability is a rather straightforward issue when only one particle is left in our model. 

Finally we mention the problem of triviality. It is a frustrating experience that the sim- 
plest models such as the 4 and QED which are used in the textbooks to demonstrate the 
renormalization of Quantum Field Theories might well be non-renormalizable if they are 
found to be trivial. In this case the study of their ultraviolet scaling behaviour serves phe- 
nomenological interest and a real ultraviolet fixed point can be achieved by asymptotically 
free models only. The one-loop ultraviolet structure of our theory turned out ot be similar 
to the conventional 4 model. This leads us to assume that our theory is not asymptotically 
free and perhaps trivial because its coupling constant which is marginal at the tree-level 
becomes irrelevant due to the one-loop contributions. This suggests the extension of the 
investigation of the antiferromagnetic vacuum to other, more involved asymptotically free 
models which may preserve their renormalizability and offer a more consistent example of a 
non-homogeneous vacuum which actually appears homogeneous in the experiments. 
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FIG. 1. The phase boundary between the paramagnetic (P), ferromagnetic (F) and the anti- 
ferromagnetic (AF) phase for C4 = 0. The two particles are degenerate along the chiral symmetric 
lines xp, Xf, Xaf- (a): The plane {m 2 LR , M\)\ (b): The plane (c2,m 2 LR ). The arrows show the 
different continuum limits at the critical point CR. 




FIG. 2. The one-loop self-energy graphs. 
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TABLES 



phase 






Z(l) 


Z(16) 


P 
F 

AF 


m LR 

-2m 2 LR 
-2m\ R - 3M 2 L 


m LR + M 2 L 
-2m 2 LR + M 2 L 
-2ml R - 2M 2 L 


1 
1 
1 


-1 + 32c 2 
-1 + 32c 2 
-1 + 32c 2 



TABLE I. The parameters of the propagator for B\ and B\§. 



phase 


K 


a 


m\ 


m|(16) 


m|(l)-m|(16) 


Pf 





-1 


> 


m R + 16fi 2 (n) 


-16// 2 (/c) 


Xp 


> 


±1 


> 


m\ 





Paf 





+1 


> 16/x 2 


m 2 R — 16fJ?(n) 


16/z 2 (/«) 


P- AF 





+1 


= 16// 2 





16/i 2 (K) 


AF 





+1 


< 16/x 2 


-2m| + 32^ 2 (K) 


16^ 2 (k) 


XAF 


> 


+1 


< 


-2m 2 R 





XF 


> 


-1 


< 


-2m\ 





F 





-1 


< 


-2m 2 R + 16^ 2 (k) 


-16/x 2 (k) 



TABLE II. The different ways of approaching the critical point of Fig. lb. Pp and Paf are the 
regions in the paramagnetic phase respectively on the left and on the right side of the chiral line 
Xp- 
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APPENDIX A: COMPUTATION OF THE EFFECTIVE POTENTIAL 

The tree-level: The tree-level effective potential for the background field 

^=$i7 1 + W, (Al) 
in the antiferromagnetic phase is the sum of the renormalized potential and the counterterms, 

^<°>(*) = \ (m 2 LR + 5ml + 3(X R + 5\)^\ 

+l(p 2 (2 d )JC(-P 2 (2 d )) + m\ R + 5ml + 3(A* + 5X)cp 2 2 M 2 2d 



2 

+ (5m 2 L + 5\lfil d )lfi 2 d^2'i + <P2*(*R + ^)($2d + 3$i$ 2 d) 

+ ^+sx { ^ i + ^ d + ml ^ (A2) 



what can be written up to a constant as 

V AFW (*) = \ (p 2 (l)/C(-P 2 (f )) + m 2 LR + 5m 2 L ) *? 

+ i (p 2 (2 rf )/C(-P 2 (2 rf )) + mi fl + 5m|) ($ 2d + ^ 



+ Afl+jA [$ 4 + ($2d + 02d)4] + 3 ( ^ + 5A)$ 2 ($2d + ^ 



The one-loop level: The next step is to obtain the one-loop contribution, 
r AF(l),^ _ 1 f d d p 



^ Fll, (*»^L /2 (^ trhi r + ^«- (p+p)2 » 

+m 2 LR + 6A fl ($i + 7 2d $ 2d )v^7 2d 
+3A J? ($ 1 + 7 2d $ 2d ) 2 ^ 



Since (7 2d ) 2 = 7 1 = 1 we can write this integral as 



(A3) 



(A4) 



/ 7 |^lndet[A(p) + P 7 2d ], (A5) 



where 
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(Mv) 
o 



A(p) + Brf 



n 

U • • 




u 




A 2\P) 




u 


n 
U 






















... o 


B\ 






... b 









... o 









... o 









and 



/ 



+ ; ; 

B 

\B 



A a (p) = G-\p) + C, 

B = 6A/J<l>i(v32d + $2<0, 

C = 3A R (2^ 2d $ 2[i + $ 2 d + $i)- 



The determinant in question is 



2 d- 



det[A(p) + 5 7 2 ] = II(^^- 52 )- 



In this manner we obtain 



V. 



AF(1) 



eff 



d d p 

2 Jp<«/2 (27r) a 



lndet[A(p) + B7 2 



1 



dp 



Y,HMp)Mp)-b 2 



2 y p <,r/2 (2vr) a=1 
In order to isolate the UV divergences it is advantageous to write 



Vffi\*) = \l 7S4E + + C) - B< 



/p< 7 r/2 (2vr) 4 ^ 



j4 2 d_1 

' /; Y, M 1 + ^ + CG * + - 52 )c a G s ] 



P<-/2 (2tt) 4 ^ 



+ 



rf 4 » 2d_1 



p<-/2 (2vr) 4 



The detailed form of flA9Q reads as 



'd^ 2 d 



p<7t/2 (2vr) d 
+m 2 l F + 3A /? (2 V 9 2d $ 2d + ^ + $ 2 ) 



(P(a)+p) 2 /C(-(P(a)+p) 2 ) 



x 
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P(a)+p) 2 JC(-(P(a)+p) 2 ) 

+™ 2 LR + SX R (2if 2d ^ 2d + $ 2 d + $ 2 ) 



(All) 



The corresponding expressions for the ferromagnetic phase are obtained by the exchange 
1 <-> 2 d , 



" e// 1 } 2 7 P < V 2 (2tt)« 

+mf R + 3A i? (2 V 9i$i + $ 2 d + $ 2 ) 

"(P(«)+p) 2 /C(-(P(a)+p) 2 ) 

+m 2 ^ + 3A i? (2 V 9 1 $ 1 + $ 2 d + $ 2 ) 

-36A 2 ,(^ 1 + $ 1 ) 2 <l> 2 d }. 

In the paramagnetic phase one sets ip a — 0, 

d d p 2 ' 



P(a)+p) a /C(-(P(a)+p) 5 



x 



(A12) 



" e// 1 ' 2 i P < V 2 (2tt)^ I 
+m| /? + 3A i? ($ 2 d + $ 2 ) 



P(a)+p) 2 /C(-(P(a)+p) 2 ) 



x 



(P(«)+p) 2 /C(-(P(a)+p) 2 ) 



+m 2 R + 3A i? ($ 2 d + $ 2 ) 



-36A|$ 2 



?*!<}■ 



(A13) 



APPENDIX B: EXPANSION AROUND THE SYMMETRICAL THEORY 

We study in this Appendix the dependence of the loop integrals on the scale parameter 
pi what controls the strength of the explicit breaking of the chiral symmetry. For this end 
we need the continuum limit for the integrals of the type 



lim/ n = lim/ d A pG™{p). 

o-»0 a~*0j p <JL. ay ' 



(Bl) 



Close to the symmetrical point the propagator has two maxima, namely in the regions 
a — 1, 16. In case of a convergent integral one expects that most of the contribution comes 
from the regions around these maxima. But due to the presence of divergencies the problem 
has to be considered more carefully. 
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The Brillouin zone B\ : The /i-dependence shows up in the antiferromagnetic phase only. 
One expands around p = 0, 

G^ n (p) = (p 2 + m 2 R {l) + baY + coV + 0(a 6 p 8 ))" (B2) 

where 

b = b + b 1 (jm) 2+K 

c = c + ci(/ia) 2+K (B3) 
are dimensionless functions. One finds 
| d 4 pG?(p) = / 



(p 2 + m 2 R {l)) r 



x 



6a 2 p 4 + ca 4 p 6 + 0(a 6 p 8 ) 
p 2 + m R (l) 



n(n + 1) fba 2 p A + ca 4 p 6 + 0(a 6 p s )\ 2 
+ 2 I p 2 + m|(l) J + 

The corrections to the usual first term are of the form 

i V+^(i))^ 1 j 

where I — k — 2m = in order to keep the dimension of each contribution. The one-loop 
integral is finite for 4 — 2(n + m) + £ — k < i.e. for n > 3. Since > the finiteness implies 
vanishing corrections and /z-independence, 

lim P d'pGUp) = l'^ d 4 p(p 2 + rh 2 R {l))-\ (B6) 

la 

for n > 3. The non- vanishing corrections arise for n — 1 and 2. 

The mass renormalization can be carried out without difficulties for fi ^ 0. In order to 
follow the renormalization of the coupling constant we need the /z-dependence for n — 2, 



-a d^p ~ r°° d A p 1 



SSy_^L (2tt)4 Gi(p) " i-oo (2tt)4 (pa + + ^ m 

where A = 2n/a and If in is a finite function of /i. Since If in depends on \i through the 
combination (/ia) 2+K and the UV divergence is logarithmic only the finite part of ( |B7|) 
becomes /i-independent in the continuum limit. 

The Brillouin zone B\§ : The mass is /i-dependent in each phase and we write rh R (16) = 
rh R (l) + rh 2 R where 

m 2 R = -16(x/i 2 (a/i) K . (B8) 
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The expansion is made around p = -P(16), 

G-^{p) = (f + rh 2 R (l) + m 2 R + bay + ca A p* + 0{aY))\ 



(B9) 



d 4 P G n M = J d 4 p 



x 



1 -n 



{p 2 +m 2 R {l)) n 

+ ba 2 p 4 + ca 4 p 6 + 0(a 6 p 8 ) 



p 2 + m 2 R (l) 



(BIO) 



n(n + 1) /fh% + ba 2 p 4 + ca 4 p 6 + 0(a 6 p 8 )\ 2 



p 2 + m|(l) 



The repetition of the argument followed in the previous case yields to the /i-independent 
finite result, 



lim/_ 2 ;^ 6 (p)=/ 



d 4 p 



oo (p 2 + m 2 R (l) + m*r 



for n > 3. For n = 2 one finds 



lim / 2 " - — ^-jG 2 Jp) 

a-W_- (2tt) 4 16W 



(27r) 4 (p 2 + m|(l)+m|) 2 + /in 



In 



A 2 



167T 2 V m|,(l) + m| 



1 ) + I fin 



In 



A 2 



16tt 2 V"m 2 R (l) 
There is a similar result for the mixed product, 



In 



^(16) 
m|(l) 



1 +/ 



/in- 



oo(27r)^p 2 +m 2 R (l))(p 2 +m 2 R (l)+m%) 



+1* 



fin 



or 



1 A A 2 m|(16) ^ m|(16) \ 

167T 2 \ m 2 ^ T7T,^.(1) / 

The regions B a , a = 2, ■ ■ ■ , 15: We can find a real number, 7, such that 

G a (p) < 7a 2 , 



r a d 4 pG n a {p) < -f n a 2n r a d 4 p. 



(Bll) 



(B12) 



(B13) 



(B14) 



(B15) 



This integral is vanishing in the continuum limit for n > 3. For n = 2 one has 
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G a (p)Gp(p) = G* a (p)G* a (p) 



1 + a/2 K a K - 2 sin 4 (p + P{a))G*Jp) 



+ sm\p + P(P))GUp + • 



(B16) 



From 



and 



fi n a n - 2 sin 4 (p + P{a)) < \ia n ~ 2 (B17) 
GUp) < la 2 (B18) 



one obtains 



fx n a n - 2 sin 4 (p + P{a))Gf{p) < /x 7 3 a n+4 (B19) 
which yields the equation 

lim J2 d'pG a (p)G p (p) = J[ a d'pGl{p)G%p) (B20) 



2a 2a 



in the continuum limit. 

So one finds the divergent part 

V 2 = V* + SV 2 , 



V 2 = P* + 5V 2 , (B21) 



with V 2 = 2V\ and 



RT) 1 In ^ (16) 

8f> 1 A m^(l)+m| ln m^(16) \ 
167T 2 V m^j(l) / ' 



(B22) 



APPENDIX C: THE FINITE PART OF THE EFFECTIVE POTENTIAL 

To compute the finite part of the effective potential we start from equation ( |A10| ), sep- 
arate the finite contributions and cancel the divergences by the counterterms. Actually the 
finite part of the potential in the continuum limit depends only the terms containing G\ et 

G\Q. 

Non-degenerate masses, (k = 0): We seek 

v£®l*) = \l ^M(P 2 + ™Ul) + C)(p 2 + m 2 R (lQ) + C)-B 2 ] 
C r d A p ( 1 1 



2 J P <A (2tt) 4 V + m|(l) p 2 + m|(16) 
C 2 r d 4 p ( 1 1 



4 J P <A (27r) 4 V(p 2 + m|j(l)) 2 (p 2 +m|(16)) 



2 7 p <a (27r) 4 (p 2 + m 2 J (l))(p 2 + m 2 ? (16))' 1 ] 
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which gives after integration 

1 



with 



($) 



1287T 2 



>nd - Kd) ln 



+ (a 2 



b 2 nd 



dnd) hi 



m|(l) 

^CL nd b nd 



d 



nd 



+d 



"* ' ' ^"2m|(l)m|(16) 
m|(16)+m|(l) ln m|(16) 



nd ~ 



"^M 1 ) - n»fl(16) m|(l) 



-2c nd (m^(16) + m^(l)) 



"6c^ — C?nd — (a n d + Kid)\J (d-nd ~ Kid) 2 + ^nd 
dnd + Kd — \J {o-nd ~ Kd) 2 + 2dnd \ 



x In 



dnd + Kd — 



-b 



nd) 



+ 2d 



nd 



V 



(C2) 



We find 



aw = rh 2 R (l) + 3A R (2v?i 6 $i6 + $f 6 + $1 

Cnd 



^16 



^2 ^F(l) 
a *l V fin 



di.V< 



<E>1 /in 



($) 

''(cD) 



$=0 



t„^: (i) (4) 



AF(1) 



FT Ft V 
a $i a $i 6 v fin 



$=0 



$=0 



m^(16) + 3A R (2y9 16 $ 16 + $ 2 6 

3A R (2y9 16 $ 16 + $ 2 6 + $ 2 ), 

72A 2 ? $ 2 ( V 9 16 + $ 16 ) 2 , 
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Degenerate masses, (/c > 0): If the two masses are degenerate then the effective potential 



is 
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In this case 
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=0 ~~ °*16 K /m 
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The effective potential and its derivatives of the ferromagnetic phase can be obtained 
from the corresponding formulae of the antiferromagnetic phase by exchanging the index 
1 < — ► 16. In the paramagnetic phase one finds 



d 2 V P{1) 

°<$- 1 V fin 



$=0 



d 2 v p{1) 



$=0 



0. 



(C8) 



Observe that the particle of the zone £>i 6 remains massless along the P-AF transition line. 
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